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We present a seareh for a narrow scalar or vector resonance decaying into Z7 with a subsequent 
Z decay into a pair of electrons or nmons. The data for this search were collected with the DO 
detector at the Fermilab Tevatron pp collider at a center of mass energy y/s = 1.96 TeV. Using 1.1 
(1.0) fb~^ of data, we observe 49 (50) candidate events in the electron (muon) channel, in good 
agreement with the standard model prediction. Prom the combination of both channels, we derive 
95% C.L. upper limits on the cross section times branching fraction (u x B) into Z7. These limits 
range from 0.19 (0.20) pb for a scalar (vector) resonance mass of 600 GeV/c^ to 2.5 (3.1) pb for a 
mass of 140 GeV/c^ 
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PACS numbers: 12.15.Ji, 13.85.Rm, 13.85.Qk, 14.70.Bh, 14.70.Hp, 14.70.Pw, 14.80.-j, 14.80.Cp 



Despite its tremendous success, the standard model 
(SM) in its current form may be a low energy approx- 
imation of a more fundamental theory. The SM does 
not describe gravity, and fundamental parameters such as 
masses and coupling constants are not derived from the 
theory. Many models exist to replace or extend the SM. 
A heavy partner of the Z boson, Z', appears in grand uni- 
fied theories, little Higgs models, models with extra spa- 
tial dimensions, and superstring theories. Scalar Higgs 
bosons, pseudo-scalar toponium, vector Z' bosons, and 
techniparticlcs could decay into the diboson final state 

This Letter presents a search for a narrow scalar or 
vector resonance decaying into Z7 using approximately 
1 fb^^ of data collected with the DO detector in col- 
lisions at -^s = 1.96 TeV at the Fermilab Tevatron col- 
lider. This analysis considers leptonic decays of the Z 
boson into electron or muon pairs. A similar search had 
been carried out by the DO collaboration using a smaller 
dataset corresponding to about 300 pb^^ 0, 

The main components of the DO detector 9] include 
a central tracking system, a sampling calorimeter, and 
a muon detection system. The central tracker consists 
of a silicon microstrip tracker (SMT) and a central fiber 
tracker (CFT), both surrounded by a 2 T superconduct- 
ing solenoid. These two subsystems provide tracking and 
vertexing for charged particles with pseudorapidities up 







FIG. 1: Diagrams for the leading-order processes which pro- 
duce candidates: (a) SM initial state radiation (ISR), (b) 
SM final state radiation (FSR) , (c) qq pair annihilation into a 
vector iy) particle which couples to the and (d) SM Higgs 
production and decay. 



to 1 77 1 K, 3 and \r]\ ~ 2.5, respectively. The pseudora- 
pidity is defined as 77 = — ln[tan(0/2)], where 9 is the 
polar angle in the DO coordinate system with the ori- 
gin at the geometrical center of the detector, the z-axis 
pointing along the beam line and the y-axis vertical up- 
ward. The liquid- argon/uranium calorimeter is spatially 
divided into the central calorimeter (CC, \ri\ < 1.1) and 
two end calorimeters (EC, 1.5 < \r]\ < 4.2). Both the 
CC and EC are segmented into electromagnetic (EM) 
and hadronic (HAD) sections. The outermost subdetec- 
tor, the muon system, consists of three layers of tracking 
detectors and scintillation trigger counters, and a 1.8 T 
toroidal iron magnet between the first two inner layers. 
The muon system is capable of providing measurements 
for particles with pseudorapidities \r]\ < 2. Arrays of 
plastic scintillators situated near the beamline in front of 
the EC are used to measure the instantaneous luminos- 
ity of the colliding proton-antiproton beams. Events in 
the electron (muon) channel are collected using a suite of 
single electron (single or dimuon) triggers, correspond- 
ing to an integrated luminosity of 1.1 (1.0) fb^^. The 
uncertainty on the luminosity is 6.1% [10|. 

A scalar model and a vector model for a particle decay 
into Z7 are used; both assume a resonance total width 
smaller than the resolution of the detector (Table In 
this narrow-width approximation, interference effects are 
negligible ll|. We use the SM Higgs boson production 
model, as implemented in pythia [12], for the scalar res- 
onance decay into the Zj final state. To model the vector 
resonance decay, we use a generic color-singlet, neutral, 
vector particle (V) implemented in madevent [l^. We 
assume Cqq is the coupling between the V and initial 
state fermions qq, with q — u 01 d and that C^-y is the 
coupling between the V and Z7, as shown in Fig. [1] A 
Z' is a good example oi a. V particle, but there is no 
model of fundamental Z' coupling to Z7, since the Z' 
has no electric charge. However, if the Z' has a compos- 
ite structure, as in technicolor models, then such a decay 
is possible. 

Electromagnetic objects such as electrons and photons 
are required to have an isolation 
and 0.15, respectively. 



of less than 0.2 
Electron candidates from the Z 
boson decay are reconstructed from EM showers in the 
calorimeter with an electron-like shower shape, that are 
required to satisfy the following criteria: have transverse 
energy Et > 15 GeV, deposit at least 90% of their energy 
in the EM calorimeter, and have tracks spatially matched 
to the EM showers. At least one electron candidate is 
required to be reconstructed in the central calorimeter, 
while the second candidate can be reconstructed either 
in the CC or EC. In addition, at least one electron must 
have Et > 25 GeV to satisfy single-electron high-i?T 
triggers. 
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Reconstruction of the Z ^ /i^ decays begins with 
a search for a pair of muons with transverse momen- 
tum pt > 15 GeV/c. To reduce the effects of the 
muon trigger px turn-on, at least one muon must have 
Pt > 20 GeV/c. Cosmic-ray background is reduced by 
rejecting muon candidates that do not originate from 
the same vertex or are reconstructed back-to-back with 
an opening angle |A0^p + A0^^ — 2tt\ less than 0.05, 
where A^^^ and A^^^ are the muon candidates sepa- 
rations in polar and azimuthal angles. The azimuthal 
angle is defined as (j) = arctan(^) so that = f points 
along the y-axis. Contamination from hadronic bb pro- 
duction is reduced by the additional requirement that 
one or both of the muon candidates are isolated from 
other activity in the calorimeter and central tracker. For 
calorimeter isolation the Et sum of calorimeter cells 
within an (77, (j)) annulus centered on the muon trajec- 
tory with R = ^(A(/))2 + (Arjy between 0.1 and 0.4 
should be less than 2.5 GeV. For central tracker isolation 
the sum of the transverse momenta of all tracks recon- 
structed within a cone centered on the muon trajectory 
with R = 0.5 should be less than 3.5 GeV/c. The trans- 
verse momentum of the muon track itself is not included 
in this sum. 

The Z"f final state is obtained by requiring an event 
to also have a photon candidate, reconstructed from an 
isolated shower in the EM part of the central calorimeter 
with the shower shape consistent with that of a photon. 
The photon candidate must be separated from both lep- 
tons from the Z boson decay by AR > 0.7. A jet with 
most of its energy carried by photons, and mis-identified 
as an electron or a photon candidate is defined as an EM- 
like jet. EM-like jet background is reduced by requiring 
that the sum of transverse momenta of all tracks, recon- 
structed within an (77, (jj) annulus centered in a photon 
candidate's trajectory, is below 1.5 GeV/c. The annulus 
has a radius R between 0.05 and 0.4. 

The three-body mass (Mu^) resolution directly affects 
the sensitivity in searching for a narrow mass resonance. 
The Mgg^ resolution is 8-18% in the muon channel and 4- 
5% in the electron channel for SM sources. To improve 
the resolution for the reconstructed three-body invari- 
ant mass in the muon channel, a Z mass constraint is 
applied to adjust the muon transverse momenta in the 
dimuon channel. A function is minimized by varying 
the muons' pt with the dimuon mass constrained to the 
Z mass; the fit has one degree of freedom. If the after 
minimization, Xmim is less than seven, the constrained fit 
is used for the Z momentum and mass. In cases where the 
constraint Xmin greater than seven, the mass constraint 
is not applied. The choice of Xmin < 7 is found by scan- 
ning values of x^ ^nd calculating the acceptance for the 
fifj.'j final state. This requirement ensures that Drell-Yan 
events will only be moved to the Z boson mass if they are 
consistent with that mass. In Monte Carlo studies this 
technique improves the Mu^ resolution substantially in 



the muon channel, as shown in Table ID 

TABLE I: Dilepton-plus-photon invariant mass resolution. 
The fourth column shows an improvement in the dimuon-plus- 
photon invariant mass resolution after the Z mass constraint 
has been applied. 

^^ee"/ ^^fi^f ^^[i^'f 

Mass (GeV/c^) resolution uncorrected corrected 

(%) resolution (%) resolution (%) 



100 


4.1 


8.3 


4.6 


200 


4.4 


8.0 


4.7 


300 


4.8 


9.9 


5.2 


400 


4.9 


11.7 


5.5 


500 


4.9 


11.8 


5.8 


600 


4.9 


13.6 


6.1 


700 


4.9 


14.6 


6.4 


800 


5.2 


16.8 


7.2 


900 


5.3 


17.6 


7.8 



To improve the analysis sensitivity in an unbiased fash- 
ion, an optimization of the photon Et and dilepton in- 
variant mass Mil selection criteria is performed with re- 
spect to S/^/S + B using a simulation. Here S is the 
number of signal events and B is the number of back- 
ground events. The simulated signal is based on MC sam- 
ples of vector and scalar resonances decaying to the Zj 
final state. The two dominant background sources, SM 
Zj and .Z-l-jets production, are considered. The back- 
ground normalization procedure is discussed in detail be- 
low. The results for optimization varied with the channel 
and resonance mass. Because the S/^/S + B has a broad 
maximum and the Mu^ distribution still had substan- 
tial discrimination, selection criteria are chosen that al- 
lowed greater acceptance. The final conditions imposed 
are£;J > 20 GeV and dilepton mass M« > 80 GeV/c^. 

The photon and lepton reconstruction efficiencies, as 
well as the total acceptances for different Z boson decay 
modes (electron and muon), depend on the spin and the 
mass of the hypothetical resonance. Electron and muon 
decay modes are also treated separately to take into ac- 
count differences in geometrical acceptance, trigger and 
reconstruction efficiencies of electrons and muons. The 
efficiency to reconstruct a pair of electrons for resonances 
with invariant masses from 120-900 GeV/c^ varies be- 
tween 60-68% (61-67%) for a vector (scalar) resonance. 
The photon reconstruction efficiency varies from 92% 
to 95% in both electron and muon channels for both 
types of resonances. The electron channel efficiencies rise 
slowly by about a factor of two up to resonance masses 
of 600 GeV/c^, whereupon they begin to drop. Simi- 
lar effects are observed in the muon channel. The muon 
identification efficiency is approximately 79% per a pair 
of muons for all resonance masses in either model. Single 
electron triggers are (99 ± 1)% efficient, and the average 
efficiency for the muon trigger requirements is (68 ± 1)%. 
Both the total efficiency of the event selection criteria 
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multiplied by the geometrical and kinematic acceptance, 
and the trigger efhciency have a noticeable mass depen- 
dence, rising from 7% to 19% for vector resonance masses 
between 120 GeV/c^ and 600 GeV/c^, and from 8% to 
20% for scalar resonance masses over the same interval in 
the electron channel. Similar effects are observed in the 
muon channel. A source of inefficiency appears at masses 
above 600 GeV/c^ where the leptons are spatially more 
coUinear and can become indistinguishable. For this rea- 
son we require the di-electron pair to be separated by 
Ai?ee > 0.6, while muon separation Ai?^^ is above 0.5. 

The two main background sources to the process un- 
der study are SM Z7 production and Z+jet production, 
where a jet is misidentified as a photon. All other back- 
ground sources were found to be negligible. To esti- 
mate the Z+}et background, we first calculate the Et- 
dependent rate, /, at which an EM-like jet is misre- 
constructed as a photon. This is done using a sample 
of events enriched with jets that satisfy the jet trigger 
requirements. The rate is the ratio of the Et spec- 
trum of photon candidates that pass all photon selec- 
tion criteria and the Et spectrum of EM objects that 
are reconstructed in the geometrical acceptance of the 
central calorimeter. The rate is further corrected for 
contamination from direct photon production. To esti- 
mate the Z+iet background two samples are used. One 
sample is our final data sample that contains events 
with a Z boson candidate and a photon candidate that 
passes all selection criteria and a sample, referred as 
to Z+FiM, consisting of events that contain a Z bo- 
son candidate plus a photon candidate that fails the 
track isolation and shower shape requirements. The lat- 
ter sample comprises of real photons and EM-like jets. 
Our final data sample contains Z + ^ events that are 
corrected by the photon identification efficiency e-y and 
Z+iet events that are corrected by /. More details on 
the calculation of the Z+jet background can be found in 
Ref. Q. The final Z+jet background is estimated to be 
4.5 ± 0.7(stat.) ± 0.6(syst.) events in the electron chan- 
nel and 4.4 ± 0.7(stat.) ± 0.6(syst.) events in the muon 
channel. The systematic uncertainty on the Z+iet back- 
ground mostly comes from the uncertainty on the photon 
efficiency and the rate at which an EM-like jet is misre- 
constructed as a photon. 

The number of SM Zj background events is estimated 
from a Zj MC sample obtained with the leading-order 
(LO) Baur event generator [l3| using CTEQ6L1 par- 
ton distribution functions (PDFs) with values of zero 
for the trilinear Zjj and ZZj couplings: NsuiZj 
Ih) = etot ■ cjsm{Zi) ■ B{Z ^ U) ■ C. Here, etot is 
the total efficiency, <7sm{Zj) ■ B{Z U) is the cross 
section times branching fraction and C is the integrated 
luminosity. We correct the LO photon Et spectrum us- 
ing the i?T-dependent ii'-factor derived from the next- 
to-leading-order (NLO) Baur event generator [l^. Ta- 
ble |TT] summarizes these numbers 0] . Using the above 



equation we obtain an estimated SM Zj contribution of 
37.4 ± 6.1(stat.) ± 2.6(syst.) events in the electron chan- 
nel and 41.6 ± 6.5(stat.) ± 2.2(syst.) events in the muon 
channel. The systematic uncertainty on the SM Zj back- 
ground comes from the uncertainty on the theoretical 
cross section, the PDFs and reconstruction efficiency. 

TABLE IL Summary of the components used to estimate 
number of SM Z'y background events. 



Parameter Electron channel Muon channel 
0.0026 ± 0.0001 0.0032 ± 0.00007 
asM (Z-f) -BiZ ^ ££) (pb) 12.85 ± 0.60 12.85 ± 0.60 
C (pb~^) 1110 ±70 1010 + 60 



The selection criteria yield 49 candidates in the elec- 
tron channel and 50 candidates in the muon channel with 
the estimated combined SM Zj plus Z+iet background 
of 41.9±6.2(stat.)±2.6(syst.) events in the electron chan- 
nel and 46.0 ± 6.6(stat.) ± 2.3(syst.) events in the muon 
channel, as shown in Tabic ITTTl 

TABLE IIP Summary of the background expectations in each 
channel and comparison with the observation. The first un- 
certainty is statistical and the second is systematic. 





Electron channel 


Muon channel 


SM Z7 


37.4 ±6.1 ±2.6 


41.6 ±6.5 ±2.2 


Z -1- jets 


4.5 ±0.7 ±0.6 


4.4 ±0.7 ±0.6 


Total background 


41.9 ±6.2 ±2.6 


46.0 ±6.6 ±2.3 


Data 


49 


50 



Figure [2] shows the distribution of the dilepton invari- 
ant mass (Mu) versus the dilepton-plus-photon invariant 
mass (M«^). The vertical band is populated by the ISR 
events where the radiated photon originates from one of 
the initial partons and the on-shell Z boson decays into 
two leptons. The Drell-Yan events cluster along the di- 
agonal band. Most of the FSR events, which would pop- 
ulate the horizontal band centered at Mu^ « Mz, are 
removed by the Mu > 80 GeV/c^ cut. 

In Fig. [21 the combined Mu^ distribution from both 
channels is compared with the SM background. Due to 
the limited available background statistics and the three- 
body mass resolution, events with > 370 GeV/c^ 
are placed into an overflow bin. Figure [4] shows the 
Mu^ distribution associated with MC signals of a vector 
particle decaying into Zj for different vector resonance 
masses. 

The observed spectrum is found to be consistent 
with SM expectations, hence limits are set on the axB for 
both vector and scalar models. The branching fraction 
for Z to ee or fijj, is accounted for in these results. 

A modified frequentist method is used to examine 
the Mu^ spectrum in the data (Fig. [3]) for discrepancies 
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FIG. 2: Invariant mass of the dilepton system vs. invariant 
mass of dilepton-plus-photon candidates. 



FIG. 4; Shape comparison of the invariant dilepton-plus- 
photon mass spectrum associated with MC signal of a vector 
particle decaying into Z7 for vector resonance masses of 120, 
180 and 260 GeV/c^ 
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FIG. 3: Invariant dilepton-plus-photon mass spectrum for 
lly data (dots), SM Z"/ background (solid line histogram) 
and Z+}et background (dashed line histogram). The shaded 
band illustrates the systematic and statistical uncertainty on 
the sum of backgrounds. 



with respect to SM sources. A Poisson log-likelihood ra- 
tio test statistic (LLR) |18] is used to compare the SM- 
only background hypothesis to one that incorporates a 
possible Zj resonance signal. The LLR incorporates sys- 
tematic uncertainties in the form of nuisance parame- 
ters that are integrated out assuming a Gaussian prior 
and a relative contribution to the signal and background 
uncertainties that is independent of the Muj invariant 
mass. When setting the limits using the LLR method 
in the combined electron and muon channels, a 2.3% re- 



construction efficiency times acceptance systematic un- 
certainty is applied to the MC signal. A 6.1% systematic 
uncertainty from luminosity and 5% PDF uncertainty are 
applied to the signals and SM Zj background. An addi- 
tional systematic uncertainty of 9% on the Z+jet back- 
ground is due to the photon efficiency and the rate at 
which an EM-like jet is misreconstructed as a photon, 
whereas an additional systematic uncertainty of 2.6% on 
the SM Z7 background is due to the theoretical cross sec- 
tion and reconstruction efficiency times acceptance. Fig- 
ures [5] and [6] show 95% C.L. exclusion curves for ax B a,s 
function of the resonance mass in the vector and scalar 
models, respectively. 

In summary, we have searched for evidence of a nar- 
row Zj resonance in eej and /i/i7 final states of pp col- 
lisions at y/s = 1.96 TeV using data collected by the 
DO detector during the 2002-2006 run of the Fermilab 
Tevatron pp Collider. We observe 49 candidate events in 
the electron channel and 50 in the muon channel, con- 
sistent with expectations from SM processes. No statis- 
tically significant evidence of a resonance decaying into 
Z"f is observed. From the combination of both channels, 
we derive 95% C.L. upper limit on cross section times 
branching fraction, which ranges from 0.19 (0.20) pb for 
a scalar (vector) resonance mass of 600 GeV/c^ to 2.5 
(3.1) pb for a mass of 140 GeV/c^. 
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FIG. 5: The observed a x B 95% confidence level limit for a 
scalar particle decaying into Z7 as a function of the scalar res- 
onance mass. The observed limit is compared to the expected 
limit for a SM Higgs decaying into Z7. The two shaded bands 
represents the 1 s.d. (dark) and 2 s.d. (light) uncertainties 
on the expected limit. 
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FIG. 6: The observed a x B 95% confidence level limit for a 
vector particle decaying into Z7 as a function of the vector res- 
onance mass. The observed limit is compared to the expected 
limit for a generic color-singlet, charge-singlet, vector particle 
decaying into Z7. The two shaded bands represents the 1 s.d. 
(dark) and 2 s.d. (light) uncertainties on the expected limit. 
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